Development and Optimization of Acriflavine-Loaded Polycaprolactone Nanoparticles Using Box–Behnken Design for Burn Wound Healing Applications

Nanoparticles are used increasingly for the treatment of different disorders, including burn wounds of the skin, due to their important role in wound healing. In this study, acriflavine-loaded poly (ε-caprolactone) nanoparticles (ACR-PCL-NPs) were prepared using a double-emulsion solvent evaporation method. All the formulations were prepared and optimized by using a Box–Behnken design. Formulations were evaluated for the effect of independent variables, i.e., poly (ε-caprolactone) (PCL) amount (X1), stirring speed of external phase (X2), and polyvinyl alcohol (PVA) concentration (X3), on the formulation-dependent variables (particle size, polydispersity index (PDI), and encapsulation efficiency) of ACR-PCL-NPs. The zeta potential, PDI, particle size, and encapsulation efficiency of optimized ACR-PCL-NPs were found to be −3.98 ± 1.58 mV, 0.270 ± 0.19, 469.2 ± 5.6 nm, and 71.9 ± 5.32%, respectively. The independent variables were found to be in excellent correlation with the dependent variables. The release of acriflavine from optimized ACR-PCL-NPs was in biphasic style with the initial burst release, followed by a slow release for up to 24 h of the in vitro study. Morphological studies of optimized ACR-PCL-NPs revealed the smooth surfaces and spherical shapes of the particles. Thermal and FTIR analyses revealed the drug–polymer compatibility of ACR-PCL-NPs. The drug-treated group showed significant re-epithelialization, as compared to the controlled group.


Introduction
Each year, acute thermal burns and injuries affect nearly half a million Americans who require medical treatment, of which approximately 40,000 need hospitalization. In the past four decades, the survival rate of burn wounds has improved amazingly; this is attributed to new treatment strategies for burn wounds, advancements and improvements in burn-care units, research on burn-wound management, and artificial skin grafting [1]. All these improvements helped us in controlling burn mortalities. Nearly 80% of burns are due to wet sources (scalds) and dry sources (flames and fire), and can be distinguished based on burn depth. A thermal injury greater than 20% of the total body surface area or greater leads to burn shock. Burn shock is characterized by fluids and protein movement from intravascular to interstitial space, increased capillary permeability, a hydrostatic pressure increase in the microvasculature, a decrease in cardiac output, and a decrease in body fluids, which leads to hypovolemia [2,3].

Codes
Factors (Independent Variables)

Z4
Primary emulsion stirring speed 3000 rpm for 3 min

Design for Optimization
Optimization is a robust approach for an accurate and precise formulation in all aspects. Design Expert (Design Expert ® version 11.0.5.0, State Ease Incorporation, Minneapolis, MN, USA) was used for the optimization of the formulation using a Box-Behnken design with three levels (low −1, medium 0, high +1) as shown in Table 2. In this study, three independent factors, including polymer concentration (X1), stirring speed (X2), and surfactant concentration (X3) was used over three dependent responses (Table 1), i.e., particle size (Y1), polydispersity index (Y2), and entrapment efficiency (Y3), for the optimization of the formulations. As shown in the Table 3, a total of 15 experimental runs were executed to acquire the optimized formulation with the statistical data, which were further processed for physicochemical characterization. The optimized formulation was selected on the basis of factor desirability over responses [22].  The prepared formulations were subjected to particle size, polydispersity, and zeta potential analyses using Zetasizer (Nano ZS, Malvern, Malvern WR14 1XZ, UK). About 1 mM NaCl was added prior to the zeta analysis to produce isotonic stability. A 100-fold dilution of the samples was performed to avoid multiple scattering using deionized water. A scattering angle of 90 • at 25 ± 1 • C for each sample was used for characterization [23].

Encapsulation Efficiency (EE)
The %EE of optimized acriflavine nanoparticle formulations were measured using slightly modified indirect methods. A specified amount of particulate suspension from each formulation was centrifuged using centrifugation tubes at 15,000 rpm for 30 min at 4 • C. The supernatant layer was collected and analyzed for free drug content at 416 nm by adding 2.5 mL of the sample in the spectrophotometer covit, and the %EE of the optimized formulation was measured using the following formula: Encapsulation efficiency (%EE) = Total acriflavine − acriflavine in supernatent Total acriflavine × 100 A standard curve was produced by preparing a series of acriflavine dilutions in deionized water. The results of all formulations were quantified using a standard curve [22].

Particles Morphology
The optimized ACR-PCL-NPs' morphology was assessed using a scanning electron microscope (SEM). Before the SEM analysis, the optimized ACR-PCL-NPs were lyophilized using a freeze-dryer (Biobase, Shandong, China). Lyophilized nanoparticles were mounted on aluminum stubs supported by adhesive tape. The SEM (Carl Zeiss Inc., Oberkochen, Germany) was operated to visualize the morphology of the ACR-PCL-NPs under high vacuum at 10 KV accelerated voltage [24,25].

Differential Scanning Calorimetry (DSC)
DSCs of pure acriflavine, PCL, a physical mixture of acriflavine and PCL, and ACR-PCL-NPs were performed to assess any change in physical state (PerkinElmer, Pyres 6.0 DSC, Waltham, MA, USA). A sample of 3 mg was used in a heating pan on a heating range of 40-300 • C and with a heating rate of 10 • C/minute, using nitrogen gas with a 20 mL/minute flow [26,27].

Attenuated Total Reflectance Fourier Transform Infrared (ATR-FTIR) Spectroscopy
ATR-FTIR was used to determine drug-interaction with components and the chemical composition of each ingredient in the formulation. Pure acriflavine, PCL, and ACR-PCL-NPs spectra were obtained using ATR-FTIR fitted with an ATR sampling cell (Thermo Scientific, Waltham, MA, USA). A resolution of 4 cm −1 was used to obtain FTIR spectra between 4000 and 400 cm −1 . Two scans of each sample were performed and spectra were obtained using OPUS 5 software.

In Vitro Release Study
An in vitro release study of optimized ACR-PCL-NPs was performed using the Franz diffusion cell method. A sample of 1 mL was applied on a cellulose membrane and it was placed between the donor and receptor compartments, and a phosphate buffer of pH 7.4 was added to the receptor compartment and was stirred at 250 rpm using a magnetic bead with the temperature maintained at 32 ± 1 • C. A sample of 1 mL was collected at different time intervals (0.5 h, 1 h, 1.5 h, 2 h, 4 h, 8 h, 12 h, 16 h, 20 h, and 24 h) and the buffer volume was replaced with 1 mL fresh dissolution media (pH 7.4). Collected samples were diluted and the acriflavine content was determined using a UV-visible spectrophotometer at 416 nm. Furthermore, kinetic release models were applied to the acriflavine release data from the optimized ACR-PCL-NPs to estimate the acriflavine release mechanism [28,29].

Histological Examination
All the experiments involving rats were conducted in accordance with NIH (USA) Care and Use of Lab Animals (NIH, 25 June 1985). The optimized formulation was tested in vivo for wound healing activity. For this purpose, 24 male rats were purchased from the National Institute of Health, weighing 250 ± 10 g. They were housed in separate cages and provided with free access to food and water for 7 days. Then, they were divided into 3 equal groups. Group 1 was a control group, group 2 was treated with an acriflavine-containing formulation, and group 3 was treated with a marketed drug (1% acriflavine). In order to anesthetize the lab animals, rats were injected intramuscularly with ketamine (40 mg/kg of body weight) and xylazine (5 mg/kg body weight). After 14 days of this treatment, rats were sacrificed to excise and collect the burn-wound sites for hematoxylin and eosin staining (H and E staining). The excised wound was washed using normal saline water and it was formalin-fixed. An approximately 3-5-µm thick section of the excised wound was stained using H and E staining to be photographed under a light microscope at 5× and 40× magnification [30]. Moreover, quantitative percent re-epithelialization was determined using the following formula: where B is a re-epithelialized skin area and C is an unclosed wound [31]. Statistical analysis was performed and the data are expressed as mean ± SD from three separate observations. For different content assays, a one way ANOVA test (p < 0.05) was used to analyze the differences among EC50. A probability of p < 0.05 was considered as significant.

Optimization of ACR-PCL-NPs
A Box-Behnken design was adopted for the ACR-PCL-NPs' optimization. The effects of the independent variables were investigated over each dependent variable using a contour plot and a 3D response surface plot. A polynomial equation was used to check the combined and individual effect of each factor on each response using Design Expert software. The quadratic effect was best applied on all factors because it offers the maximum effect both individually and in combination. The Design software managed the analysis of variance (ANOVA) of individual responses, and the results indicate the model-fitting for the data sets, as shown in Table 4 (linear, 2FI, quadratic). The factors such as polymer concentration (X1), stirring speed (X2), and surfactant concentration (X3) were evaluated at three levels (low −1, medium 0, high +1) to formulate ACR-PCL-NPs.
The selection of optimized ACR-PCL-NPs was made on the basis of specific criteria, i.e., the small size of the nanoparticles, by Design Expert.
A total of 15 formulations with three central points were fed into the Design Expert software. After analysis, the particle size (Y1) of formulation F1 was found to be 302.4 nm, which was the smallest of all 15 formulations, as shown in Table 3. Similarly, the particle size of formulation F14 was 780.4 nm, the highest of all formulations. F1 had the lowest, at 0.202% PDI (Y2), while F13 has the highest, at 0.349% PDI. The encapsulation efficiency (Y3) of formulation F2 was the lowest, i.e., 67.7%, while formulation F15 was the highest, at 81.8%, as shown in Table 3. The coefficient of correlation (R 2 ) values were in the range of 0.8699 to 0.9864, which shows the highest coherence of fit data with high value PRESS, as described in Table 4.

Particle Size, Zeta Potential, and PDI
The hydrodynamic particle size of ACR-PCL-NPs was examined. The obtained values of the ACR-PCL-NPs were the average of three independent measurements. The particle size and PDI of the optimized ACR-PCL-NPs were 469.2 ± 5.6 nm and 0.270 ± 0.19, respectively. A small particle size means a large surface area for release, and a low PDI ensures the homogeneity of the particles [22]. Zeta potential is an analytical technique used to measure the surface charge of NPs in colloidal dispersions. The zeta potential magnitude provides signs of colloidal stability. NPs with zeta potential values between >+25 and <−25 mV have high degrees of stability; low zeta potential values may lead to coagulation, flocculation, and aggregation due to van der Waals forces [32,33]. The zeta potential of ACR-PCL-NPs were found between −3.75 and −7.96 mV, which were considered to be near zero ( Figure 1). The results demonstrate no considerable variation in the zeta potentials of the formulations, which is associated with the uncharged chemical nature of PCL, as already reported in the literature [34,35].

Particle Size, Zeta Potential, and PDI
The hydrodynamic particle size of ACR-PCL-NPs was examined. The obtained values of the ACR-PCL-NPs were the average of three independent measurements. The particle size and PDI of the optimized ACR-PCL-NPs were 469.2 ± 5.6 nm and 0.270 ± 0.19, respectively. A small particle size means a large surface area for release, and a low PDI ensures the homogeneity of the particles [22]. Zeta potential is an analytical technique used to measure the surface charge of NPs in colloidal dispersions. The zeta potential magnitude provides signs of colloidal stability. NPs with zeta potential values between >+25 and <−25 mV have high degrees of stability; low zeta potential values may lead to coagulation, flocculation, and aggregation due to van der Waals forces [32,33]. The zeta potential of ACR-PCL-NPs were found between −3.75 and −7.96 mV, which were considered to be near zero ( Figure 1). The results demonstrate no considerable variation in the zeta potentials of the formulations, which is associated with the uncharged chemical nature of PCL, as already reported in the literature [34,35]. Contour plots and 3D surface plots were used to check the effect of independent variables on the ACR-PCL-NPs, as shown in Figure 2. The PCL concentration (X1) had a positive effect on the particle size of nanoparticles. An increase in PCL concentration, from 1 g to 2 g (F1 and F10), led to an increase in particle size, from 302.4 nm to 730 nm of ACR-PCL-NPs, as shown in the Table 3. Our results are in accordance with Lepeltier and coworkers' study, where they state that an increase in polymer concentration causes an increase in the nanoparticles' size [36] [37]. The stirring speed of the secondary emulsion (X2) had a negative impact on the size of the prepared nanoparticles. As we increased the stirring speed of the external phase from 8000 rpm (F11) to 12,000 rpm (F5), a decrease in particle size from 670 nm to 499 nm was observed, which was similar to what was demonstrated by Ibraheem et. al. [20]. This decrease in particle size may be due to the provision of high levels of energy to the particulate system, as high stirring speeds cause the dismantling of large droplets of the second emulsion into small droplets [38]. As we increased the stirring speed of the formulations, a gradual decrease in the particle size of Contour plots and 3D surface plots were used to check the effect of independent variables on the ACR-PCL-NPs, as shown in Figure 2. The PCL concentration (X1) had a positive effect on the particle size of nanoparticles. An increase in PCL concentration, from 1 g to 2 g (F1 and F10), led to an increase in particle size, from 302.4 nm to 730 nm of ACR-PCL-NPs, as shown in the Table 3. Our results are in accordance with Lepeltier and co-workers' study, where they state that an increase in polymer concentration causes an increase in the nanoparticles' size [36,37]. The stirring speed of the secondary emulsion (X2) had a negative impact on the size of the prepared nanoparticles. As we increased the stirring speed of the external phase from 8000 rpm (F11) to 12,000 rpm (F5), a decrease in particle size from 670 nm to 499 nm was observed, which was similar to what was demonstrated by Ibraheem et al. [20]. This decrease in particle size may be due to the provision of high levels of energy to the particulate system, as high stirring speeds cause the dismantling of large droplets of the second emulsion into small droplets [38]. As we increased the stirring speed of the formulations, a gradual decrease in the particle size of the formulations was detected. PVA should be used in optimum concentrations; low concentrations of PVA (X3) lead to coalescence, while high concentrations of PVA can cause clumping of the nanoparticles [39]. The combined effect of both polymer concentration (X1) and surfactant concentration (X3) was positive for the particle size, and the effect of stirring speed in the external phase (X2) was negative on the particle size. the formulations was detected. PVA should be used in optimum concentrations; low concentrations of PVA (X3) lead to coalescence, while high concentrations of PVA can cause clumping of the nanoparticles [39]. The combined effect of both polymer concentration (X1) and surfactant concentration (X3) was positive for the particle size, and the effect of stirring speed in the external phase (X2) was negative on the particle size.
The polynomial equation (Equation 1) indicates that X1 and X3 have positive effects on the particle size (Y1) of the ACR-PCL-NPs (p < 0.0001). The model's F-value, 16.28, and p-value, 0.0034, suggests the model's significance. The lowest p-value and highest F-value indicates the best fit of the model. The Prob F-value was 0.005. All the independent variables showed a significant effect on the ACR-PCL-NPs. The R 2 and adjusted R 2 values were in close range, which shows the model's robust credibility.
The polynomial equation (Equation (1)) indicates that X1 and X3 have positive effects on the particle size (Y1) of the ACR-PCL-NPs (p < 0.0001). The model's F-value, 16.28, and p-value, 0.0034, suggests the model's significance. The lowest p-value and highest F-value indicates the best fit of the model. The Prob F-value was 0.005. All the independent variables showed a significant effect on the ACR-PCL-NPs. The R 2 and adjusted R 2 values were in close range, which shows the model's robust credibility.

Effect of Independent Variables on PDI (Y2)
The polymer concentration (X1) has shown a high positive impact on the PDI of the formulations. The higher the PCL concentration that was used, the higher the PDI that was observed, as shown in Table 3. The stirring speed of the external phase (X2) had a slightly negative impact on the PDI, and this was in accordance with Scholz et al.'s findings [40]-that an increase in stirring speed can decrease both the particle size and PDI of the system. Polymer concentration was the prominent factor in determining the formulation of the PDI [41]. The surfactant concentration (X3) had a very slightly negative impact on the PDI of a formulation, as shown in Table 3 and Figure 3. F13, with a PVA concentration of 0.3%, and F14, with a PVA concentration of 0.7%, had 0.399 and 0.363 PDI, respectively, which is in agreement with the work of Tefas et al. [42], where they stated that PVA has a negative effect on the PDI of the formulation when used as a surfactant. The F-value of the model was 40.28, which made the model significant. The p-value 0.0004 indicates only a 0.041% chance of noise to occur in the model.

Effect of Independent Variables on PDI (Y2)
The polymer concentration (X1) has shown a high positive impact on the PDI of the formulations. The higher the PCL concentration that was used, the higher the PDI that was observed, as shown in Table 3. The stirring speed of the external phase (X2) had a slightly negative impact on the PDI, and this was in accordance with Scholz et al.'s findings [40]-that an increase in stirring speed can decrease both the particle size and PDI of the system. Polymer concentration was the prominent factor in determining the formulation of the PDI [41]. The surfactant concentration (X3) had a very slightly negative impact on the PDI of a formulation, as shown in Table 3 and Figure 3. F13, with a PVA concentration of 0.3%, and F14, with a PVA concentration of 0.7%, had 0.399 and 0.363 PDI, respectively, which is in agreement with the work of Tefas et al. [42], where they stated that PVA has a negative effect on the PDI of the formulation when used as a surfactant. The Fvalue of the model was 40.28, which made the model significant. The p-value 0.0004 indicates only a 0.041% chance of noise to occur in the model.  Polydispersity Index (Y2): The higher F-value (40.28) and the lower p-value (0.0004) suggest that the model is significant and stable. The R 2 value and the adjusted R 2 value were 0.9864 and 0.9619, respectively. The model significance is a sign of a model with low noise.

Effect of Independent Variables on Entrapment Efficiency (%EE) (Y3)
The effect of different factors over %EE (Y3) can be seen in the contour plot and 3D surface plot, as given in Figure 4. The polynomial equation (Equation (4)) showed the positive impact of polymer concentration (X1) and surfactant concentration (X3) on %EE, as demonstrated by Ibraheem et al. [20,36]. Formulations F2 and F13, with %EE of 67.7 and 78.9, respectively, have polymer concentrations of 1g and 2 g, respectively. An increase in polymer concentration had a direct influence on %EE, i.e., increasing the % EE [38]. Similarly, the surfactant concentration (X3) had a positive influence on the % EE. An increase in PVA concentration (X3) will lead to an increase in %EE, because PVA has the tendency to encapsulate the drug in the matrix [43]. The stirring speed (X2) showed a slightly negative impact on the %EE. An increase in stirring speed (X2) led to a decrease in %EE, which is accordance with Bozena et al.'s work, where they stated that increases in stirring speed will decrease the drug content of particles [44]. Similarly, JK Patel's study also showed a decrease of %EE with an increase of stirring speed, which is an accordance with our results [45].
The higher F-value (40.28) and the lower p-value (0.0004) suggest that the model is significant and stable. The R 2 value and the adjusted R 2 value were 0.9864 and 0.9619, respectively. The model significance is a sign of a model with low noise.

Effect of Independent Variables on Entrapment Efficiency (%EE) (Y3)
The effect of different factors over %EE (Y3) can be seen in the contour plot and 3D surface plot, as given in Figure 4. The polynomial equation (Equation 4) showed the positive impact of polymer concentration (X1) and surfactant concentration (X3) on %EE, as demonstrated by Ibraheem et al. [20,36]. Formulations F2 and F13, with %EE of 67.7 and 78.9, respectively, have polymer concentrations of 1g and 2 g, respectively. An increase in polymer concentration had a direct influence on %EE, i.e., increasing the % EE [38]. Similarly, the surfactant concentration (X3) had a positive influence on the % EE. An increase in PVA concentration (X3) will lead to an increase in %EE, because PVA has the tendency to encapsulate the drug in the matrix [43]. The stirring speed (X2) showed a slightly negative impact on the %EE. An increase in stirring speed (X2) led to a decrease in %EE, which is accordance with Bozena et al.'s work, where they stated that increases in stirring speed will decrease the drug content of particles [44]. Similarly, JK Patel's study also showed a decrease of %EE with an increase of stirring speed, which is an accordance with our results [45].
The polynomial equation (Equation (5)) suggests the positive effect of both PCL (X1) and PVA (X3) on %EE, and the negative impact of stirring speed (X2) on %EE (p < 0.001). The model's F-value, 24.06, suggests the model is significant. The high F-value and low p-value (0.0013) delineate the fitness and reliability of the model. The R 2 and adjusted R 2 value were 0.9774 and 0.9368, respectively, which are close to each other. The significance and low noise in the model enhances the credibility of the model.

Preparation of Optimized ACR-PCL-NPs
Optimized ACR-PCL-NPs were prepared using same method as followed for all 15 trial formulations. The optimized formulation was formulated on the basis of values given by the Design Expert software after the analysis of 15 trial formulations for independent variables. Optimization was based on the small particle size, low polydispersity index, and high entrapment efficiency. After the analysis of the trial formulations, a formulation with a desirability of 1 was selected as the optimized formulation, having all the above-discussed properties. The software provided actual values for the independent and predicted values of the dependent variables of the optimized formulation. The actual values calculated after optimization for the PCL concentration (X1), the stirring speed of the external phase (X2), and the PVA concentration (X3) were 1.5g, 12,000 rpm, and 0.5%, respectively. The predicted values of the responses for optimized formulation were 451.9 nm for Y1, 0.253 for Y2, and 72.34% for Y3. When the optimized formulation of ACR-PCL-NPs was prepared using actual values of independent variables, the responses' actual values were 469.2 ± 5.6 nm for Y1, 0.270 ± 0.19 for Y2, and 71.9 ± 5.32% for Y3. The responses' actual values were close to the predicted values. The optimized ACR-PCL-NPs were processed for characterization, for an in vitro release study, and for a kinetic model-fitting study.

Surface Morphology
A scanning electron microscope (SEM) was used for the morphological analysis of the ACR-PCL-NPs prepared by using the double-emulsion solvent evaporation method. The surface texture, shape, inter particulate bridging, and smoothness of the ACR-PCL-NPs were evaluated. The SEM images show that the ACR-PCL-NPs were smooth on the surface with a spherical shape. The smoothness of the ACR-PCL-NPs' surface supports the assumption that encapsulated active drug-release may be due to matrix erosion [46]. Additionally, bridging between particles can be seen in SEM images ( Figure 5), and these bridges may be associated with the sticky nature of the PVA used in these formulations. Being sticky in nature, it is very difficult to completely remove the PVA, even after particle washing [47]. 3.2.6. Differential Scanning Calorimetry Figure 6 represents the DSC thermograms of pure acriflavine, PCL, physical mixtures of PCL and acriflavine, and the ACR-PCL-NPs. The DSC thermogram of pure acriflavine represents a sharp endothermic peak at 180 °C, which shows the crystalline nature of the acriflavine. The PCL DSC thermogram showed a sharp peak at 61.2 °C. The physical mixture of acriflavine and PCL showed two distinctive peaks, one at 61.2 °C for PCL, and another at 180 °C for acriflavine. However, the peak of acriflavine in the ARR-PCL-NPs was absent because of the encapsulation of acriflavine in the PCL matrix in an amorphous form [5,48]. A single peak, at 61.2 °C, was observed in the ACR-PCL-NPs' thermogram. The melting temperature (Tm) of PCL and ACR-PCL-NPs remained the same (61.2 °C), as shown in Figure 6. The degree of crystallinity (Xc) of the ACR-PCL-NPs decreased (PCL = 93% and ACR-PCL-NPs = 75.5%). The enthalpy of fusion (Δ Hm) of PCL and ACR-PCL-NPs were 126.5 J/g and 103 J/g, respectively. The decrease in Xc and ΔHm in the ACR-PCL-NPs was due to a lack of homogenous distribution and structure uniformity [49]. 3.2.6. Differential Scanning Calorimetry Figure 6 represents the DSC thermograms of pure acriflavine, PCL, physical mixtures of PCL and acriflavine, and the ACR-PCL-NPs. The DSC thermogram of pure acriflavine represents a sharp endothermic peak at 180 • C, which shows the crystalline nature of the acriflavine. The PCL DSC thermogram showed a sharp peak at 61.2 • C. The physical mixture of acriflavine and PCL showed two distinctive peaks, one at 61.2 • C for PCL, and another at 180 • C for acriflavine. However, the peak of acriflavine in the ARR-PCL-NPs was absent because of the encapsulation of acriflavine in the PCL matrix in an amorphous form [5,48]. A single peak, at 61.2 • C, was observed in the ACR-PCL-NPs' thermogram. The melting temperature (Tm) of PCL and ACR-PCL-NPs remained the same (61.2 • C), as shown in Figure 6. The degree of crystallinity (Xc) of the ACR-PCL-NPs decreased (PCL = 93% and ACR-PCL-NPs = 75.5%). The enthalpy of fusion (∆ Hm) of PCL and ACR-PCL-NPs were 126.5 J/g and 103 J/g, respectively. The decrease in Xc and ∆Hm in the ACR-PCL-NPs was due to a lack of homogenous distribution and structure uniformity [49].  Figure 7 represents the FTIR spectra of acriflavine, PCL, and the ACR-PCL-NPs. Th characteristic peaks of acriflavine were found at 3207 cm −1 , which indicates the stretchin of symmetric C-NH, at 1324 cm −1 , which indicates the asymmetric C-N stretching vibra tions, at 1130 cm −1 , which was due to C-H stretching, at 1383 cm −1 , which was for charac teristic cyclic C-N-C vibrations, and at 1171 cm −1 , which was for the acriflavine CH group. The CH-CH3 group can be confirmed by the peak at 968 cm −1 ; the 929.81-cm −1 pea indicates CH3 vibrations. PCL is an aliphatic polymer that shows characteristic peaks be tween 2860-3000 cm −1 due to C-H stretching. The 1719-cm −1 peak was due to carbony group C=O stretching vibrations. The 1044-cm −1 -1298-cm −1 peaks were due to saturated ester C-O stretching. The 735-cm −1 peak represents C-H bending vibrations. The peaks o PCL and ACR-PCL-NPs were similar in structure; only the decrease in the peak of th ACR-PCL-NPs, from 1723 cm −1 to 1721 cm −1 , was due to C=O stretching of unsaturated ester. The peaks of the ACR-PCL-NPs were not clearly observed in the formulation FTIR spectra. This might be due to the presence of acriflavine in molecular dispersion in a pol ymer matrix [5]. The FTIR spectra do not show any potential interactions or chemical in compatibilities between acriflavine and PCL in the nanoparticles.  Figure 7 represents the FTIR spectra of acriflavine, PCL, and the ACR-PCL-NPs. The characteristic peaks of acriflavine were found at 3207 cm −1 , which indicates the stretching of symmetric C-NH, at 1324 cm −1 , which indicates the asymmetric C-N stretching vibrations, at 1130 cm −1 , which was due to C-H stretching, at 1383 cm −1 , which was for characteristic cyclic C-N-C vibrations, and at 1171 cm −1 , which was for the acriflavine CH 3 group. The CH-CH 3 group can be confirmed by the peak at 968 cm −1 ; the 929.81-cm −1 peak indicates CH 3 vibrations. PCL is an aliphatic polymer that shows characteristic peaks between 2860-3000 cm −1 due to C-H stretching. The 1719-cm −1 peak was due to carbonyl group C=O stretching vibrations. The 1044-cm −1 -1298-cm −1 peaks were due to saturated-ester C-O stretching. The 735-cm −1 peak represents C-H bending vibrations. The peaks of PCL and ACR-PCL-NPs were similar in structure; only the decrease in the peak of the ACR-PCL-NPs, from 1723 cm −1 to 1721 cm −1 , was due to C=O stretching of unsaturated ester. The peaks of the ACR-PCL-NPs were not clearly observed in the formulation FTIR spectra. This might be due to the presence of acriflavine in molecular dispersion in a polymer matrix [5]. The FTIR spectra do not show any potential interactions or chemical incompatibilities between acriflavine and PCL in the nanoparticles.

In Vitro Release Study
The solution of acriflavine showed a 96.5 ± 2.5% release in 1 h, while the optimized ACR-PCL-NPs showed 85.43 ± 5.38% in 24 h. The ACR-PCL-NPs showed a burst release followed by a sustained release, revealing a biphasic release profile (Figure 8). The first burst release was due to the surface drug on ACR-PCL-NPs, and the controlled release was due to the diffusion of acriflavine from the NPs into the release medium [50]. The acriflavine diffusion from the ACR-PCL-NPs was due to its low molecular weight and high solubility in the dissolution medium [51]. The core drug diffuses much more slowly, due to a longer diffusion path. The sustained release was favorable for a prolonged therapeutic effect. The release data of the optimized ACR-PCL-NPs were incorporated into different release models and evaluated for release mechanisms, as shown in Figure 9. To understand the process of drug release from nanoparticles, the data obtained were analyzed using four in vitro release models: zero-order, first-order, the Higuchi model, and the Korsmeyer-Peppas model.

In Vitro Release Study
The solution of acriflavine showed a 96.5 ± 2.5% release in 1 h, while the optimized ACR-PCL-NPs showed 85.43 ± 5.38% in 24 h. The ACR-PCL-NPs showed a burst release followed by a sustained release, revealing a biphasic release profile (Figure 8). The first burst release was due to the surface drug on ACR-PCL-NPs, and the controlled release was due to the diffusion of acriflavine from the NPs into the release medium [50]. The acriflavine diffusion from the ACR-PCL-NPs was due to its low molecular weight and high solubility in the dissolution medium [51]. The core drug diffuses much more slowly, due to a longer diffusion path. The sustained release was favorable for a prolonged therapeutic effect.

In Vitro Release Study
The solution of acriflavine showed a 96.5 ± 2.5% release in 1 h, while the optimized ACR-PCL-NPs showed 85.43 ± 5.38% in 24 h. The ACR-PCL-NPs showed a burst release followed by a sustained release, revealing a biphasic release profile (Figure 8). The first burst release was due to the surface drug on ACR-PCL-NPs, and the controlled release was due to the diffusion of acriflavine from the NPs into the release medium [50]. The acriflavine diffusion from the ACR-PCL-NPs was due to its low molecular weight and high solubility in the dissolution medium [51]. The core drug diffuses much more slowly, due to a longer diffusion path. The sustained release was favorable for a prolonged therapeutic effect. The release data of the optimized ACR-PCL-NPs were incorporated into different release models and evaluated for release mechanisms, as shown in Figure 9. To understand the process of drug release from nanoparticles, the data obtained were analyzed using four in vitro release models: zero-order, first-order, the Higuchi model, and the Korsmeyer-Peppas model. The release data of the optimized ACR-PCL-NPs were incorporated into different release models and evaluated for release mechanisms, as shown in Figure 9. To understand the process of drug release from nanoparticles, the data obtained were analyzed using four in vitro release models: zero-order, first-order, the Higuchi model, and the Korsmeyer-Peppas model. In 1961, Higuchi proposed the Higuchi model to describe drug-release from matrix systems. The equation is: where Qt is the amount of drug released from the united surface area "A" in time "t". C is the initial drug concentration, Cs is the drug solubility in the matrix medium, and D is the drug diffusivity in a matrix [52]. The Higuchi model describes the drug release from a matrix as a diffusion process based on Fick's law: when release kinetics are according to this model, the drug release from the particles will mainly be controlled by a diffusion process in a polymer matrix. To dig further into the mechanism of drug-release kinetics, the data obtained were analyzed using the Korsmeyer-Peppas model; the particles are considered spherical according to the morphological data [53]. The Korsmeyer-Peppas model is a semi-empirical model used for the analysis of different dosage forms' release data, and the equation is given as follows: where Mt is the amount of drug at a time "t", M∞ is the amount of drug at infinite time, K is the constant, and "n" is the release exponent.
To determine the n value in the Korsmeyer-Peppas model, only the Mt/M∞ ≤ 0.6 release fraction will be used. If the value of n is 0.43, it indicates that the dosage form releases the drug following Fickian diffusion. If the value is between 0.43 and 0.85, it will indicate the release of the drug by diffusion as well as by polymer chain relaxation [54]. In our case, release kinetics fitted the Higuchi model with an R 2 value of 0.9917, which indicates that the active release of the drug from the dosage form was mainly controlled by a diffusion process. The data fitted in the Korsmeyer-Peppas model showed an R 2 value of 0.8497 with an n value of 0.5285, which indicates drug release by two phenomena: drug diffusion and polymer chain relaxation. The Higuchi model is the perfect model to describe the drug dissolution profile of modified-release formulations, such as matrix In 1961, Higuchi proposed the Higuchi model to describe drug-release from matrix systems. The equation is: where Qt is the amount of drug released from the united surface area "A" in time "t". C is the initial drug concentration, Cs is the drug solubility in the matrix medium, and D is the drug diffusivity in a matrix [52]. The Higuchi model describes the drug release from a matrix as a diffusion process based on Fick's law: when release kinetics are according to this model, the drug release from the particles will mainly be controlled by a diffusion process in a polymer matrix. To dig further into the mechanism of drug-release kinetics, the data obtained were analyzed using the Korsmeyer-Peppas model; the particles are considered spherical according to the morphological data [53]. The Korsmeyer-Peppas model is a semi-empirical model used for the analysis of different dosage forms' release data, and the equation is given as follows: where Mt is the amount of drug at a time "t", M∞ is the amount of drug at infinite time, K is the constant, and "n" is the release exponent.
To determine the n value in the Korsmeyer-Peppas model, only the Mt/M∞ ≤ 0.6 release fraction will be used. If the value of n is 0.43, it indicates that the dosage form releases the drug following Fickian diffusion. If the value is between 0.43 and 0.85, it will indicate the release of the drug by diffusion as well as by polymer chain relaxation [54]. In our case, release kinetics fitted the Higuchi model with an R 2 value of 0.9917, which indicates that the active release of the drug from the dosage form was mainly controlled by a diffusion process. The data fitted in the Korsmeyer-Peppas model showed an R 2 value of 0.8497 with an n value of 0.5285, which indicates drug release by two phenomena: drug diffusion and polymer chain relaxation. The Higuchi model is the perfect model to describe the drug dissolution profile of modified-release formulations, such as matrix tablets and transdermal systems with water-soluble drugs. The release-and kinetic-model studies were in agreement with the study of Miladi et al. [50].

Wound Re-Epithelialization
The percentage of the total epithelial covering formation on the wounded area at a specific time point is called % wound re-epithelialization. On day 1, the re-epithelialization of the wound was 0%. The rate of re-epithelialization on specific time points was greater in the formulation-treated group, as compared to the marketed drug-treated group and the control group. On day 3, the formulation drug-treated group had 17.14% ± 3.68% re-epithelialization, and the marketed drug-treated group and the control group had 10.12 ± 3.9% and 4.13% ± 2.11%, respectively (p < 0.001). On day 7, a significant difference was observed: the formulation drug-treated group was leading with 38.12% ± 8.15%, and the marketed drug-treated group and the control group had 31.10% ± 10.32% and 14.22% ± 15.54% re-epithelialization of a wound, respectively (p < 0.05). On day 14, the formulation drug-treated group achieved an almost-complete wound re-epithelialization (98.12% ± 0.32%), the marketed drug-treated group had 84.21% ± 3.12% (Table 5), and the control group had 65.22% ± 16.07% wound re-epithelialization (p < 0.001). Table 5. Rate of re-epithelialization of wounds at the 3rd, 7th, and 14th day (n = 6).

Tissue Growth and Gross Histology
The H-and E-stained slides of the controlled group, the marketed drug-treated group, and the formulation-treated group of rats were prepared for histological evaluation. The H-and E-stained slides were visualized and photographed using a light microscope at 5× and 40× magnification, which allowed us to visualize the recovered wound area and the formation of epidermis and dermis. The formulation-treated group showed an increase in burn-wound thickness, as compared to the marketed drug-treated and control groups. The formulation drug-treated group completed epidermal regrowth, with intact epidermises on the 14th day, while the marketed drug-treated group still had acute inflammatory infiltrations. The gap between the epidermis and dermis was small in the formulationtreated and marketed drug-treated groups, as compared to the control group. Blood vessels were smaller and higher in number in the formulation-treated group as compared to the marketed drug-treated group and the control group. The control group was much further behind in the regeneration process of the dermis and epidermis, and the persisting redness of the burn wound indicated inflammation and irregular granulation ( Figure 10). The H and E staining study was also quantified by using a total skin-thickness measurement on a universal testing machine, showing an increase in the thickness of the formulation-treated group (1849 ± 240 µm, p < 0.01) and the marketed drug-treated group (1682 ± 340 µm, p < 0.05), compared to the untreated group (1152 ± 150µm). Hair follicles and rete pegs were not observed in the central wound, even though some were detected on the wound edges, which was likely due to the regeneration and migration of endogenous cells.

Conclusions
In this study, the developed optimized ACR-PCL-NPs showed desirable physicochemical properties. Nanoparticles were prepared using a double-emulsion solvent evaporation method, using PCL polymer for the encapsulation of the drug and polyvinyl alcohol as a stabilizing agent. A Box-Behnken design was applied to 15 formulations using Design Expert 11 ® . The developed ACR-PCL-NPs showed a particle size less than 780 nm, and a drug-entrapment efficiency of about 69.1%, with a spherical shape and a nearly neutral surface charge. The in vitro release was found to be biphasic, initially featuring a burst release for 1 h, followed by a sustained release for 24 h. The %age re-epithelization of the formulation drug-treated group (98.12 ± 0.32) was higher, as compared to the marketed drug-treated group (84.21 ± 3.12) and the untreated group (65.22 ± 16.07). Therefore, from all the above results, it can be concluded that optimized ACR-PCL-NPs could be a productive delivery system for acriflavine in burn-wound healing.

Conclusions
In this study, the developed optimized ACR-PCL-NPs showed desirable physicochemical properties. Nanoparticles were prepared using a double-emulsion solvent evaporation method, using PCL polymer for the encapsulation of the drug and polyvinyl alcohol as a stabilizing agent. A Box-Behnken design was applied to 15 formulations using Design Expert 11 ® . The developed ACR-PCL-NPs showed a particle size less than 780 nm, and a drug-entrapment efficiency of about 69.1%, with a spherical shape and a nearly neutral surface charge. The in vitro release was found to be biphasic, initially featuring a burst release for 1 h, followed by a sustained release for 24 h. The %age re-epithelization of the formulation drug-treated group (98.12 ± 0.32) was higher, as compared to the marketed drug-treated group (84.21 ± 3.12) and the untreated group (65.22 ± 16.07). Therefore, from all the above results, it can be concluded that optimized ACR-PCL-NPs could be a productive delivery system for acriflavine in burn-wound healing.